arXiv:1503.02770vl [nucl-th] 10 Mar 2015 


Electromagnetic deuteron form factors in point 
form of relativistic quantum mechanics 


N. A. Khokhlov 

Komsomolsk-na-Amure State Technical University 

Abstract 

A study of the electromagnetic structure of the deuteron in the framework 
of relativistic quantum mechanics is presented. The deuteron form factors de¬ 
pendencies on the transferred 4-momentum Q up to 7.5 Fm” 1 are calculated. 
We compare results obtained by different realistic deuteron wave functions stem¬ 
ming from Nijmegen-I, Nijmegen-II, JISP16, CD-Bonn, Paris and Moscow (with 
forbidden states) potentials. The nucleon form factors parametrization consis¬ 
tent with the modern experimental analysis was used as the input data. 
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1 Introduction 

In the Born approximation of a one-photon exchange mechanism the elastic ed scat¬ 
tering observables are directly expressed through the electromagnetic (EM) deuteron 
form factors (FFs) [TH3]- Therefore this process allows to extract the EM FF depen¬ 
dencies on the transferred 4-momentum Q in the spacelike region. Relativistic effects 
may be essential even at low Q 013 ]. There are different relativistic models of the 
deuteron EM FFs HHEj- 

We apply point-form (PF) relativistic quantum mechanics (RQM) to treat the 
elastic electron-deuteron scattering in a Poincare-invariant way. Concepts of the RQM 
and an exhaustive bibliography are presented in the review by Keister and Polyzou |9|. 
The PF is one of the three forms of RQM proposed by Dirac [10]. The other two are 
the instant form and front form. These forms are associated with different subgroups 
of the Poincare group which may be free of interactions. General method for putting 
interactions in generators of the Poincare group was derived in m- It was shown that 
all the forms are unitary equivalent [12] . Though each form has certain advantages, 
there are important simplifying features of the PF m In the PF all generators of the 
homogeneous Lorentz group (space-time rotations) are free of interactions. Therefore 
only in the PF the spectator approximation (SA) preserves its spectator character in 
any reference frame [T4lH6| . For an electromagnetic NN process it means that the 
NN interaction does not affect the photon-nucleon interaction, and therefore the sum 
of the one-particle EM current operators is the sum of the one-particle EM current 
operators in any reference frame (r. f.). Two equivalent SAs of EM current operator 
for composite systems in PF RQM are derived in Refs. [KHITo] . The PF SA was 
applied to calculate EM FFs of deuteron, pion, nucleon 0 HMD with reasonable 
results. 

Present paper is an extension of our previous investigations where we described 
elastic NN scattering up to 3 GeV of laboratory energy [23], bremsstrahlung in the 
pp scattering pp —> ppj [23], photodisintegration of deuteron y.D —► np [2TH2(!1 and 
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exclusive electrodisintegration of deuteron m- Here we show that the developed 
approach is applicable to the elastic eD scattering. 


2 Potential model in PF of RQM 


In PF of RQM a system of two particles is described by the wave function, which is 
an eigenfunction of the mass operator M. We may represent this wave function as 
a product of the external and internal parts. The internal wave function |%) is also 
an eigenfunction of the mass operator and for system of two nucleons with masses 
m i ss m 2 ~ m = 2 m\m 2 /(jn\ + m 2 ) satisfies the equation 


M|X) = 


2\/q 2 + m 2 + V int 


lx) = M|x), 


(1) 


where Vi n t is an operator commuting with the full angular momentum operator and 
acting only through internal variables (spins and relative momentum), q is a momen¬ 
tum operator of one of the particles in the center of mass frame (relative momentum). 
Rearrangement of m gives 

[q 2 +mV] |x) = <? 2 |x), ( 2 ) 


where V acts like V)nt only through internal variables and 


q 


2 



(3) 


Eq. © is identical in form to the Schrodinger equation. The only relativistic cor¬ 
rection here is in the deuteron binding energy that must be changed by the effec¬ 
tive value 2.2233 MeV instead of the experimental 2.2246 MeV. It is easy to show. 
Let £ be the deuteron binding energy. Then M = 2m — e and for deuteron state 

q 2 = -m 2 = —me (l — -r~)- Comparing with the nonrelativistic relationship 

q 2 = —me we can identify factor (l — -A -) as the essential relativistic correction. 
There is no similar correction in the scattering region because q 2 = mE\ a b/2 is the 
precise relativistic relationship ( Ei a b is the laboratory energy), that is used in the 
partial wave analysis. The change is negligible for our problem. Therefore this formal 
identity allows us to use non relativistic deuteron wave functions in our calculations. 


3 eD elastic scattering 

We give only results of PF RQM necessary for our calculation, in notation of Ref. [15] . 
We use formalism of m for calculation of the matrix elements of the EM current 
operator. 

We consider the pn system and neglect difference of neutron and proton masses. 
Let pi be the 4-momentum of nucleon i, P = (P , P) = pi + P 2 be the system 4- 
momentum, M be the system mass and G = P/M be the system 4-velocity. The 
wave function of two particles with 4-momentum P is expressed through a tensor 
product of external and internal parts 

|P,x) = E/i 2 |P)<SS|x), (4) 

where the internal wave function \x) satisfies Eqs. tU-©. The unitary operator 

2 

U12 = Pi 2 (G,q) = pP[s i ;a(p,:/ m ) _1 «(G)a:((3'i/m)] (5) 

»=1 

is the operator from the ’’internal” Hilbert space to the Hilbert representation space 
of two-particle states m- P[s; u] is the SU(2) representation operator corresponding 
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to the element u SSU(2). s are generators of the representation. In our case of 
spin s = 1/2 particles, we deal with the fundamental representation, i.e. s; = -icr, 
(er = ((T x ,(Ty,cr~) are the Pauli matrices) and D[s;u] = u. The momenta of the 
particles in their c.m. frame are 


Qi = L[a(G)} x pi, (6) 

where L[a(G)\ is the Lorentz transformation to the frame moving with 4-velocity G. 
Matrix a(g) = ( g° + 1 + a • g)/y / 2(g° + 1) corresponds to a 4-velocity g. 

The external part of the wave function is defined as 

(G\P') = ^-G°S 3 (G-G'). (7) 


Its scalar product is 

(P"\P') = J f^-(P"\ G )(G\P') = 2 V / M' 2 + P' 2 (5 3 ( P" - P'), (8) 

where G°(G) = y/l + G 2 . The internal part |x) is characterized by momentum 
q = qi = —q 2 of one of the particles in the c.m. frame. 

According to the Bakamjian- Thomas procedure m interaction appears in 4- 
momentum P = GM, where M is sum of the free mass operator Mf ree and of 
the interaction Vi nt : M = Mf ree + V) nt (Eq. {TJ). The interaction operator acts 
only through internal variables. Operators V)nt and V (and therefore M and Mf ree ) 
commute with spin operator S (full angular momentum) and with 4-velocity operator 
G. Generators of space-time rotations are free of interaction. Most non-relativistic 
scattering theory formal results are valid for case of two particles [S]. For example in 
the c.m. frame the relative orbital angular momentum and spins are coupled together 
as in the non-relativistic case. 

The deuteron wave function \Pi,Xi) is normalized as follows 

(P f ,Xf\Pi,Xi) = 2P°5 3 (P i - P/)<X/|Xi). (9) 

There is a convenient r. f. for calculation of current operator matrix elements m 
(it coincides for elastic ed scattering with the Breit r. f.). For all EM reactions with 
two nucleons this r. f. is defined by condition: 


G f + Gj — 0, (10) 

where G/ = P f/Mo, G, = P,;/M o are final and initial 4-velocities of the deuteron 
and Mjj is its mass. The matrix element of the current operator is US]: 

(P f , X/1 ( x )\PiiX%) = 2(M / M l ) 1/2 exp(*(P / - P i )x){ X f\''f{HXi ), (H) 


where j M (h) defines action of current operator in the internal space of the NN system. 


2 (MjMf) 1 / 2 k 

(il di + A^fy) 2 2Md 


( 12 ) 


is vector-parameter [T5] (0 < h < 1), k is momentum of photon in r. f. (flT)l) . Mi = 
Mf = Md are masses of initial and final NN system (deuteron). 

The internal wave function of deuteron is 


I Xi) = \ Uj ( r )l Z ’ X J = IM J )r 

1 = 0,2 


( 13 ) 
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with normalisation (Xi|Xi) = 1- We use the momentum space wave function: 

I Xi) = - u M)\1^'AMj) q, 

^ 1= 0,2 


(14) 


where 


u(q) = u 0 (q) 



dr sin (qr)u(r), 


w(q) = u 2 {q) = 




sin (qr) 


— cos(gr) 
qr 


w(r). 


(15) 

(16) 


Transformations from the Breit r. f. m to the initial c. nr. frame of the NN 
system and to the final one are boosts along vector h (axis z). Projection of the total 
deuteron angular momentum onto the z axis does not change for these boosts. The 
initial deuteron in the Breit r. f. moves in direction opposite to the h. Its internal 
wave function with spirality Aj is 


|A») = - Y, Mq)\l, 1; 1, Mj = —A*). (17) 

5 1= 0,2 

Wave function of the final deuteron with spirality A f is 

|A/) = - E u l (q)\l,l-,l,M J = A f ), (18) 

^ i=0,2 


Usual parametrisation of the EM CO matrix element for the deuteron (spin 1 
particle) is [21121 [28]: 


(AP®Pf) 1 / 2 (P f ,x f \J^\Pi, Xi) 

= - {ci(o 2 Kc; ■«.) - G,w 2 ) (g ' A ^‘' AP) } <p» + pjf) 

- G 2 (Q 2 )l(?((} • A P) - (}"((, • AP)1, (19) 

where ( a-b) = a°b°— (a-b), form factors Gi(Q 2 ), i = 1,2, 3 are function of Q 2 = —AP 2 , 
A P = P f -Pi. 

In the Breit r. f. P/ = -Pi, P° = Pf = P° = M D /V 1 - h 2 , A P = (0,2P/), 
P? + P$ = (2 P°, 0 ), P f /P° = h, P f = hM D /y/l=h?, A P 2 = -Ah 2 M 2 /(1 - ft 2 ), 
Q 2 = -AP 2 , h 2 = (h-h). 


(Xf\j°N\ X i) = -Gr(Q 2 )(r • 0 + 2G 3 (Q 2 ) ^ f \ h ^' h) 

+ G 2 (Q 2 )[$(e r h)-$*(&-h)], (20) 

(X/U(h)|x<) = G 2 (Q 2 )fo(€} ■ h) - CMi ■ h)] = G 2 (Q 2 )[h x fc x £}]]. (21) 

It can be shown H5!> that these expressions are equivalent to choosing j v as: 


j°(h) = Gc(Q 2 ) + ~JY~J^Gq(Q 2 ) 
j(h) = -^=L=G M (Q 2 )(h x J), 



(h-J) 2 


( 22 ) 

(23) 


where J is total angular momentum (spin) of the deuteron. Gc, Gq, Gm are its 
charge monopole, charge quadruple and magnetic dipole FFs. 
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Spiral polarizations of the deuteron in the initial and final states are 

= 


^/ A = 


(o,±i,-*,o)/V5 

(-<2/2,0,0, Po)/M d = (-h, 0,0, l)/vr^7? 
(0, =Fl, -z, 0)/V2 

(<2/2, o, o, p 0 )/m d = (h, o, o, i)/ 


(A = ±) 

(A = 0), 

(A = ±) 

(A = 0). 


Polarization of the virtual photon is 

A = f (0, =Fl, —i, 0)/\/2 

1 ( 1 , 0 , 0 , 0 ) 

FFs Gi are expressed as 

Gq = G\ + t^vGq, 

g q = g 1 -g m + (i + v)g 3 , 

2 h 2 


(A = ±) 

(A = 0). 


(24) 

(25) 

(26) 


Gi = G c - 


3(1 — h 2 ) 


Gq, 


(27) 


G, = Go 1 - 


— Gc(l — h 2 ) + Gm(1 — h 2 ), 


where = <2 2 /4M|, = h 2 /(l — h 2 ). 

At Q 2 = 0 we have Gc = Gi, Gq = G\ — Gm + G 3 . Form factors Gq( 0) = e, 
Gm(0) = [idc/ZMd and Gq = Qdc/Mq give charge, magnetic and quadruple mo¬ 
menta of deuteron. 

Denoting helicity amplitudes as j\ f \. = (A/| (e* • j M (h)) |Aj), we arrive at: 


Joo (Q) = Gc + r 


4 h 2 


;Gq, 


j+-(Q 2 ) = j° + (<2 2 ) = G c - 

j+o(<2 2 ) + j(j"_(<2 2 ) 


.3 1 - h 2 

2 h 2 


3 1 - h 2 
h 


Gq, 


VT^h? 


G m 


and 


jUQ 2 ) = j-o(Q 2 )« hUQ 2 ) = Jo+(Q 2 )- 

The deuteron FFs are associated with unpolarised structure functions 

A(<2 2 ) = Gc{Q 2 ) + g 7 7G 2 / (<2 2 ), 

B(Q 2 ) = lv(i + v)G 2 M (Q 2 ). 


(28) 

(29) 

(30) 

(31) 

(32) 

(33) 


These quantities are extracted from the elastic eD scattering with unpolarised par¬ 
ticles. A tensor polarisation observable f 2 o(<2 2 ,<0 is usually used as an additional 
quantity needed for definition of all three FFs. 

In the present paper we use SA of the EM CO of m without expanding it in 
powers of h and we calculate matrix elements in the momentum space. Therefore 
calculating (13T1) we use a following expansion of / M (h) « jg A (h) [27j 


JU h) = (1 + (A 2 ■ s 2 )) (£f + (Cj* • 8l ))Ii(h) 

+ (1 + (Ar • sO) (B% + (C£ • s 2 )) I 2 (h), (34) 

where A,;, B 1 /', C/ are some vector functions of h and q (q,9,(f>). In the spherical 
coordinate system (q, 9, <j>) dependence of these functions on <j> appears as e ±lm ^ (to = 
0,1, 2). The </> is analytically integrated giving trivial equalities in (1311) . 
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4 Results 

In our calculation we use as an input the momentum space deuteron wave functions 
and nucleon EM FFs. The momentum space deuteron wave functions stemming from 
Nijmegen-I (Nijml), Nijmegen-I (Nijmll) [30], JISP16 [31], CD-Bonn [32], Paris [33] , 
Argonnel8 [34] (momentum space deuteron wave function is from [35) ) and Moscow 
(with forbidden states) [32] potentials are shown in Figs. [H We use two versions 
of Moscow type potential: Moscow06 [21] and Moscowl4. Last one was derived by 
author in the same manner as in )22i but with deuteron asymptotic constants fitted to 
describe static deuteron form factors. Parameters of both Moscow potentials may be 
requested from the author (e-mail: nikolakhokhlov@yandex.ru). For all but JISP16 
the S -wave functions change sign at q ~ 2 Fm' 1 , and D -wave functions change sign 
at q ss 6 — 8 Fm^ 1 . The S- and D -wave functions of Argonnel8, Paris and Nijmll 
are close for q < 5 Fm -1 . The 5-wave functions of CD-Bonn and Nijml are close 
for q < 5 Fm' 1 . Wave functions of the JISP16 decrease rapidly after about 2 Fm -1 
without changing sign. Our results for deuteron EM FFs are presented in Table [T] 
and in Figs. [2] ED Q] Results for Argonnel8, Paris and Nijmll are close reflecting 
closeness of their wave functions for q < 5 Fm -1 . Nijml and CD-Bonn give more 
distinct results. Our calculations show that Gm changes sign for all potentials at 
rather low Q that is not seen experimentally. Nevertheless CD-Bonn and Nijml give 
good results for Gm for Q < 7 Fm. Moscow potentials give the best description 
of charge form factor Gc ■ The essential factor that influences our calculations is 


Table 1: Static deuteron form factors. Two values through slash are relativistic 
calculation/nonrelativistic calculation. 



Gm( 0) = 

Gq(0) = M^Qd 

Exp 

1.7148 

25.83 

Nijml 

1.697/1.695 

24.8/24.6 

Nijmll 

1.700/1.695 

24.7/24.5 

Paris 

1.696/1.694 

25.6/25.2 

CD-Bonn 

1.708/1.704 

24.8/24.4 

Argonnel8 

1.696/1.694 

24.7/24.4 

JISP16 

1.720/1.714 

26.3/26.1 

Moscow06 

1.711/1.699 

24.5/24.2 

Moscowl4 

1.716/1.700 

26.0/25.8 


the nucleon FF dependencies on the momentum transferred to the individual nucleon 
Qp ~ Qn 7^ Q 2 - These FFs are measured at discrete values of Qf n while we 
need dependencies on Qi. In our calculations we take phenomenological nucleon FF 
dependencies on Qf from [53]. It should be noted that neutron EM FFs are extracted 
from experimental data of 2 H(e, e'n)p and other processes with deuteron and triton 
using various models of these possesses and nuclei. Therefore these FFs are model 
dependent. We see good general correspondence of the theory and experiment for 
Q < 5 Fm” 1 . Discrepancies for larger Q are comparable with differences of results for 
different potentials. Model calculations [54] show that meson exchange currents may 
give significant effect in EM processes with ?zp-system. We do not take into account 
these currents. However it is not clear how these currents may be agreed with short 
range part of the NN interaction of quantum chromodynamic origin. Besides the EM 
FFs of nucleons are not described by meson degrees of freedom at intermediate and 
high energies [55] , 

To complete this line of our investigation, we plan to calculate neutron EM FFs 
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compatible with Moscow potential model which has not been used for the extraction 
of these FFs. 
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U/(q), fm 



Figure 1: Momentum space deuteron wave functions used in the calculations. Same 
legend for S- and D -wave functions. 
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Figure 2: Deuteron form factor Gc as a function of Q. Data are from compilation 
of [3] calculated from A, B and t 2 o data of [36U52] . 



Figure 3: Deuteron form factor Gq as a function of Q. Same legend as Fig. [2] 
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Figure 4: Deuteron form factors Gm as a function of Q. Same legend as Fig. [2] 





10 7 
10 8 



